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Abstract 

This note contributes to the point calculus of persistent ho- 
mology by extending Alexander duality to real- valued func- 
tions. Given a perfect Morse function / : §™ +1 — > [0, 1] and 
a decomposition = U U V such that M = U n V is an 
n-manifold, we prove elementary relationships between the 
persistence diagrams of / restricted to U, to V, and to M. 
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1 Introduction 

Persistent homology is a recent extension of the classical the- 
ory of homology; see e.g. 0. Given a real-valued func- 
tion on a topological space, it measures the importance of a 
homology class by monitoring when the class appears and 
when it disappears in the increasing sequence of sublevel 
sets. A technical requirement is that the function be tame, 
which means it has only finitely many homological criti- 
cal values, and each sublevel set has finite rank homology 
groups. Pairing up the births and deaths, and drawing each 
pair of values as a dot (a point in the plane), we get a multiset 
which we refer to as the persistence diagram of the function. 
It is a combinatorial summary of the homological informa- 
tion contained in the sequence of sublevel sets. If we substi- 
tute reduced for standard homology groups, we get a slightly 
modified reduced persistence diagram of the function. 

As between homology groups, we can observe relation- 
ships between persistence diagrams. A prime example is 
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Poincare duality, which says that the p-th and the (n — p)- 
th homology groups of an n-manifold are isomorphic. More 
precisely, this is true if homology is defined for field coeffi- 
cients, which is what we assume throughout this paper. The 
extension to functions says that the diagram is symmetric 
with respect to reflection across the vertical axis; see Q. 
Here, we change the homological dimension of a dot from p 
to n — p whenever we reflect it across the axis. This paper 
contributes new relationships by extending Alexander dual- 
ity from spaces to functions. To state our results, we assume 
a perfect Morse function, / : § n+1 — > [0, 1], which for the 
sphere has no critical points other than a minimum and a 
maximum, and a decomposition of the (n + 1) -dimensional 
sphere into two subsets, § n+1 = U U V, whose intersection 
is an n-manifold, M = U n V. Our first result says that the 
reduced persistence diagrams of / restricted to U and to V 
are reflections of each other. We call this the Land and Water 
Theorem. Our second result relates land with shore. Ignor- 
ing some modifications, it says that the persistence diagram 
of / restricted to M is the disjoint union of the diagram of 
/ restricted to U and of its reflection. The modifications be- 
come unnecessary if we assume that the minimum and max- 
imum of / both belong to a common component of V. We 
call this the Euclidean Shore Theorem. 

In the example that justifies the title of this paper, and the 
names of our theorems, we let U be the planet Earth, not 
including the water and the air. To a coarse approximation, 
U is homeomorphic to a 3-ball, sitting inside the Universe, 
which we model as a 3-sphere, S. The function we consider 
is the negative gravitational potential of the Earth, which is 
defined on the entire Universe. The sea is then a sublevel 
set of this function restricted to V, which is the closure of 
S — U. With these definitions, our results relate the persis- 
tence diagram of the gravitational potential restricted to the 
Earth with the shape of the sea as its water level rises. Also, 
the Euclidean Shore Theorem applies, unmodified, expand- 
ing the relationship to include the sea floor, which is swept 
out by the shoreline as the water level rises. 



Besides developing the mathematical theory of persistent 
homology, there are pragmatic reasons for our interest in the 
extension of Alexander duality to functions. Persistence has 
fast algorithms, so that the bulk of the work is often in the 
construction of the space and the function for which we com- 
pute persistence. A point in case is the analysis of the bio- 
logical process of cell segregation started in |6|. Modeling 
the process as a subset of space-time, the function of interest 
is time which, after compactifying space-time to S 4 , has no 
critical points other than a minimum and a maximum. The 
subset U of § 4 is a union of cells times time, whose boundary 
is a 3-manifold. We can represent U by a 1-parameter family 
of alpha complexes, whose disjoint union has the same ho- 
motopy type; see e.g. J5] Chapter III]. However, the bound- 
ary of that disjoint union is not necessarily a 3-manifold. Us- 
ing our Euclidean Shore Theorem, we can compute the per- 
sistence diagram of the function on the 3-manifold without 
ever constructing the 3-manifold. 

Outline. Sections [2] and [3] introduce the necessary back- 
ground on homology and persistent homology. Sections [4] 
and|5]present our two results. Section|6]concludes the paper. 

2 Homology 

Starting with a brief introduction of classical homology 
groups, we present the relevant background on Alexander 
duality and Mayer- Vietoris sequences. More comprehensive 
discussions of these topics can be found in textbooks of al- 
gebraic topology, such as Q H) ■ 

Background. The p-dimensional homology of a topolog- 
ical space, X, is a mathematical language to define, count, 
and reason about the p-dimensional connectivity of X. There 
are different but essentially equivalent theories depending on 
the choices one makes in the representation of the space, the 
selection of cycles, and the meaning of addition. For our 
purpose, the most elementary of these theories will suffice: 
a simplicial complex, K, that triangulates X, formal sums of 
p-simplices with zero boundary as p-cycles, and adding with 
coefficients in a field, F. Most algorithms on homology as- 
sume this model, in particular the ones developed within per- 
sistent homology, which is defined only for field coefficients. 
In this model, we call a formal sum of p-simplices a p-chain, 
a p-cycle if its boundary is empty, and a p-boundary if it 
is the boundary of a (p + 1) -chain. The p-boundaries form 
a subgroup of the p-cycles, which form a subgroup of the 
p-chains: B p C Z p C C p . The p-th homology group is the 
quotient of the p-cycles over the p-boundaries: H p = Z p /B p . 
Its elements are sets of p-cycles that differ from each other 
by p-boundaries. More fully, we denote the p-th homology 
group by H p (/v, F), or by H p (X, F) to emphasize that the 



group is independent of the simplicial complex we choose to 
triangulate the space. However, we will fix an arbitrary field 
F and drop it from the notation. For field coefficients, H p (X) 
is necessarily a vector space, which is fully described by its 
rank, /3 p (X) = rankH p (X) such that H p (X) ~ F^ x >. This 
rank is called the p-th Betti number of X. Finally, we will 
often drop the dimension from the notation by introducing 
the direct sum, H(X) = pgZ H p (X). 

Besides standard homology, we will frequently use re- 
duced homology groups, H p (X), which are isomorphic to 
the non-reduced groups except possibly for dimensions p = 
0, —1. To explain the difference, we note that /?o(X) counts 
the components of X. In contrast, /?o(X) = rankflo(X) 
counts the gaps between components or, equivalently, the 
edges that are needed to merge the components into one. 
Hence, /?o(X) = (3o(X) — 1, except when X is empty, in 
which case /?o(X) = Aj(X) = 0. To distinguish this case 
from a single component, we have (X) equal to 1 if X = 

0, and equal to 0, otherwise. Furthermore, we use relative ho- 
mology, which is defined for pairs of spaces, Y C X. Taking 
a pair relaxes the requirement of a chain to be called a cy- 
cle, namely whenever its boundary is contained in Y, which 
includes the case when the boundary is empty. We write 
H p (X, Y) for the p-th relative homology group of the pair, 
and /? p (X, Y) = rankH p (X, Y) for the p-th relative Betti 
number. As before, we will suppress the dimension from the 
notation by introducing H(X, Y) = © peZ H p (X, Y). 

As examples, consider the (n + 1) -dimensional sphere, 
§ = §" +1 , and a closed hemisphere, iCS, which is a ball 
of dimension n + 1. In standard homology, we have /3o (§) = 
(3q(M) = f3 n+ i(§) = 1 while all other Betti numbers are 
zero. In reduced homology, we have /3„ + i(§) = 1 while all 
other reduced Betti numbers are zero. In particular, /? P (H) = 
for all p. In relative homology, we have /3 n+ i(§, H) = 1 
while all other relative Betti numbers are zero. In particular, 
/3o(S,EI) = 0, which may be confusing at first but makes 
sense because every point on the sphere can be connected by 
a path to a point in the hemisphere and is thus a 0-boundary. 

Alexander duality. Recall that a perfect Morse function is 
one whose number of critical points equals the sum of Betti 
numbers of the space. For a sphere, this number is 2: a min- 
imum and a maximum. Throughout the remainder of this 
section, we assume a perfect Morse function /:§—>• [0, 1], 
whose values at the minimum and the maximum are and 

1. We also assume complementary subsets whose union is 
the sphere, U U V = §, and whose intersection is an n- 
manifold, U D V = M. We assume that the restriction of 
/ to M is tame, which implies that its restrictions to D and 
V are also tame. For each i 6 I, we write St = / _1 [0,i] 
for the sublevel set of /, and U t = U n S t , V t = V n § t , 
Mt =MnS; for the sublevel sets of the restrictions of / to 
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U, V, M. Similarly, we write S* = / -1 [t, 1] for the super- 
level set, and U* = U n 8*, V* = V n S*, M* = M n S* for 
the superlevel sets of the restrictions. 

The basic version of Alexander duality is a statement 
about two complementary subsets of the sphere; see J5] page 
424]. More specifically, it states that H g (U) and H P (V) are 
isomorphic, where q = n — p, except for p = and q = 
when the O-dimensional group has an extra generator. This 
implies 

A,(V) = /3„(U) + 1, (1) 
/3 p (V) = /3 g (U), (2) 
/3„(V) = /3 (U)-1, (3) 

for 1 < p < n — 1. We will also need the version that deals 
with two complementary subsets of the (n + 1) -dimensional 
baU, B n+1 ; see J9] page 426]. Let t G (0, 1) denote a reg- 
ular value of /|m- Note first that § t is homeomorphic to 
jjn+i gy excision, the homology groups of Ut relative to 
Ut fl f~ 1 (t) are isomorphic to those of U relative to U*. 
Alexander duality states that H 9 (U,U') and H p (V t ) are iso- 
morphic, where q+p = n, as before, except forp = 0, when 
Vt has an extra component. This implies 

A)(V t ) = /?„(U,U*) + 1, (4) 
/? P (V t ) - /yU,U 4 ), (5) 

for 1 < p < n. 

Mayer-Vietoris. We can connect the homology groups of 
U and V with those of M and § using the Mayer-Vietoris se- 
quence of the decomposition. This sequence is exact, mean- 
ing the image of every map equals the kernel of the next 
map. Counting the trivial homology groups, the sequence is 
infinitely long, with three terms per dimension: 

. . . -> H p+ i(S) ->■ H P (M) H P (U) © H P (V) 

->H P (S) ->..., 

where the maps between homology groups of the same 
dimension are induced by the inclusions. The only non- 
trivial homology groups of § are in dimensions and n + 1, 
with ranks /3 (§) = /3 n +i(§) = 1- It follows that for 
1 < p < n — 1, the groups defined by M, U, V are sur- 
rounded by trivial groups. This implies that the groups of M 
and the direct sums of the groups of U and V are isomorphic. 
For p = and p = n, the non-trivial homology groups of § 
prevent this isomorphism, and we get 

A)(M) = /3 (U) + A)(V)-1, (6) 
/3 p (M) = /yU) + /3p(V), (7) 
/3JM) = /3„(U) +/?„(¥) + !, (8) 



for 1 < p < n — 1. We also have a Mayer-Vietoris sequence 
for the spaces § t , U t , V f and M t : 

. . . -> H p+ i(St) -> Hp(Mt) -> Hp(Ut) © Hp(V t ) 

->• Hp(St) -> . . . . 

For < t < 1, §t is a point or a closed (n + 1) -dimensional 
ball, and its only non-trivial homology group is of dimension 
0, with rank /3 (§t) = 1. It follows that for 1 < p < n, 
the groups defined by Mt, Ut, Vt are surrounded by trivial 
groups, which implies that the groups of M t and the direct 
sums of the groups of Ut and ¥ t are isomorphic. We thus get 

A)(Mt) = /8 (U t ) + A)(V t )-l, (9) 
^(M*) = /3 p (Ut) + /3p(Vt), (10) 

for 1 < p < n. We may also consider the relative homology 
groups, again connected by a Mayer-Vietoris sequence: 

...-> H p+ i(S,S') H p (M,M t ) 
-> Hp(U, U*) © Hp(V, V 4 ) ->• Hp(S, S*) -> . . . . 

For < t < 1, the only non-trivial homology group of the 
pair (S, §') is in dimension n + 1, which implies 

/3 p (M,M 4 ) = ^(U.U'H/^V.V), (ID 
/3„(M,M*) = /3„(U,U*) + / 8„(V,V*) + 1, (12) 

for < p < n — 1. 

Separating manifold. Combining Alexander duality with 
Mayer-Vietoris, we get relations between the homology of 
the sublevel sets of M and U. More specifically, we get 

/3p(M) = /yU)+/3,(U), (13) 
f3 p (M t ) = /yUt)+/3,(U,U*), (14) 
/3p(M, M f ) = ^(U,U*)+^(U t ), (15) 

for all < p < n. Here, we combine ©, (fTJ, (O with ((T), 
©, © to get (fOJ. Similarly, we combine ©, (JTOj with ©, 
(O to get ( fT4l i. Finally, we exploit the symmetry between U 
and V and combine dTTJ, OS with ©, © to get <[T3J. Note 
that (O and O imply /3p(M t ) = /3 g (M, M f ), which is a 
consequence of Lefschetz duality; see J9)- 

An example. We illustrate the above relationships with an 
example. Let M be a 2-dimensional torus in § 3 . Accord- 
ingly, U and V are the two solid tori that decompose the 
3-sphere and intersect in M. This is sketched in Figure [TJ 
where we assume that U is the part of space surrounded by 
M, while V is the space surrounding M. The only non-zero 
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Betti numbers of the solid torus are /?q (U) = (3i (U) = 1 . 
We thus get 

/3 (M) = /3 (U)+/3 2 (U) = 1, 
^(M) - ^(UH/MU) = 2, 
/3 2 (M) = /3 2 (U)+/3 (U) + 1 = 1 

from ( f]~3T >- These are the correct Betti numbers of the 2- 
dimensional torus. Next, choose t so that the sublevel set of 
U is half the solid torus. Its only non-zero Betti number is 




Figure 1 : The level set denned by t splits the solid torus into two 
halves. 

/?o(Ut) = 1, and the only non-zero relative Betti number of 
the pair is /?i (U, U* ) = 1. We thus get 

/3 (M t ) = /3 (U t )+/3 2 (U,U*) = 1, 
/3 x (M t ) = ^(UO + ZM^U') = l, 
(3 2 (M t ) = /3 2 (U t )+A,(U,U 4 ) = 

from (fl4l . These are the correct Betti numbers of the sub- 
level set of M. Finally, we get 

/3 (M, M*) = /9 2 (U t )+A>(U,U t ) = 0, 
^(M,M*) = ^(UO+ACU.U*) = 1, 
/3 2 (M,M') = /3 (U t )+/? 2 (U,U') = 1 

from $15[ . These are the correct relative Betti numbers of the 
pair (M,M*). 

3 Persistence 

Starting with a brief introduction of persistent homology, we 
review its recent combinatorial expression as a point calcu- 
lus; see 15] and 0][2) f° r more comprehensive discussions. 

Background. We take the step from homology to persis- 
tent homology by replacing a space with the sequence of 
(closed and open) sublevel sets of a function on the space. 



To explain this, let X be compact and g : X — > [0, 1] tame. 
As before, we write X t = g^ 1 [0, t] for the sublevel set and 
X* = 1] for the superlevel set defined by t. Since 

g is tame, it has finitely many homological critical values, 
So to s m . We assume w.l.o.g. that s = and s, n = 1, 
else, we add them as additional critical values. We inter- 
leave the Si with homological regular values U, such that 
so < to < si < . . . < t m -i < s m . Taking the homology 
of closed and open sublevel sets at the regular values, we get 
the filtration 

Xo — > ■ ■ ■ — > X m _i — > X m — > X m+ i X 2m , 

where X, is H(X tj ), if < i < m - 1, H(X), \fi = m, and 
H(X, X* 2 " 1 ^), if m < i < 2m. For notational convenience, 
we add the trivial groups, X_i = and X2 TO +i = 0, at the 
beginning and end of the filtration. The maps connecting the 
homology groups are induced by the inclusions X s C X t and 
X* C X s , whenever s < t. The maps compose and we write 
gi j : Xi — > Xj. Reading the filtration from left to right, we 
see homology classes appear and disappear. To understand 
these events, we say a class a S X^ is born at X^ if it does not 
belong to the image of gj-i,,. If furthermore gij (ai) belongs 
to the image of gj-ij but gij-i(a) does not belong to the 
image of gi-i j-i, then we say a dies entering Xj. Since the 
filtration starts and ends with trivial groups, every homology 
class has well-defined birth and death values. Every event 
is associated with the immediately preceding homological 
critical value, namely with Si if the event is at X,;, and with 
s m -i if the event is at X m+ i + i, for < i < m — 1. Events 
at X m are associated with s m = 1. We represent a class by 
a dot in the plane whose two coordinates mark its birth and 
its death. More precisely, the coordinates signal the increase 
and decrease of Betti numbers, and the dot represents an en- 
tire coset of classes that are born and die with a. All dots 
have coordinates in [0, 1], by construction. 

Note that we use each Si twice, once during the first pass 
in which the sublevel set grows from to X, and then again 
during the second pass in which X minus the superlevel set 
shrinks back from X to 0. When we collect the dots, we 
separate the passes by drawing each coordinate axis twice. 
The result is a multiset of dots, which we refer to as the per- 
sistence diagram of the function, and denote as Dgm(g), or 
Dgm p (g) if we restrict ourselves to the dots representing p- 
dimensional classes; see Figure [2] where we label each dot 
with the dimension of the classes it represents. We further 
distinguish four regions within the diagram: the ordinary, 
the horizontal, the vertical, and the relative subdiagrams, de- 
noted as Ord(g), Hor(g), Vcl(g), and Rel(g). For example, a 
dot belongs to the ordinary subdiagram if its birth and death 
both happen during the first pass; see again Figure|2] 

Reduced persistence diagrams. It is sometimes conve- 
nient to use reduced instead of standard homology groups 
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Figure 2: The persistence diagram of the function restricted to the 
solid torus (dots drawn as circles), and to the complementary solid 
torus (dots drawn as squares), as sketched in Figure Q] The white 
dots belong to the reduced diagram, the black dots belong to the 
standard diagram, and the shaded dots belong to both. 

in the filtration. There are small differences caused by the 
(—1) -dimensional class, which exists if the space is empty. 
As a first step, we introduce reduced versions of the rela- 
tive homology groups, which are isomorphic to the standard 
relative groups, for all p. To define them, let w be a new 
(dummy) vertex, write ui ■ X* for the cone of ui over X*, 
and let X^ = X U u> ■ X* be the result of gluing the cone 
along X* to X. Then the reduced relative homology group of 
the pair is H P (X,X') = H p ( X D- For example, if X* = 0, 
then uj forms a separate component, so that the reduced Betti 
number is equal to the number of components of X, just as 
/?o(X, X*). The reduced persistence diagram is now defined 
as before, but for the filtration of reduced homology groups. 
To describe this, we suppress the dimension and write X^ 
for H(X ti ), if < i < m, for H(X), if i = m, and for 
H(xL 2m ~ l ), if m < i < 2m. The resulting sequence of re- 
duced homology groups, from Xo to )<2 m , is connected from 
left to right by maps induced by inclusion. Finally, we de- 
fine Dgm(t/) by matching up the births and the deaths and 
by drawing each pair as a dot in the plane. Similar to before, 
we write Dgm p (g) when we restrict ourselves to homology 
classes of dimension p. 

For dimension p > 1, the reduced diagrams are the same 
as the standard ones, simply because the groups are the same. 
More formally, the persistence diagrams are the same be- 
cause the two filtrations form a commutative diagram with 
vertical isomorphisms: 

Xg -> X? -> ... -> X p 2m -> 

X^ -> Xl -+ ... -> XL -> 0; 

see the Persistence Equivalence Theorem J5] page 159]. This 
is not true for p = — 1, where the standard diagram is empty, 
while the reduced diagram contains a single dot marking the 
transition from an empty to a non-empty space. To describe 
the difference for p = 0, we call a dot (it, w) in the 0-th 



standard persistence diagram extreme if no other dot has a 
first coordinate smaller than u and a second coordinate larger 
than it). Here, we assume for simplicity that no two dots in 
the 0-th diagram share the same first coordinate or the same 
second coordinate. 

We note that only dots in the horizontal subdiagram can be 
extreme. The non-extreme dots also belong to the reduced 
diagram of g, while the extreme dots exchange their coordi- 
nates to form new dots in the reduced diagram; see Figure 
[2] The way the coordinates are exchanged will be important 
later, so we describe the details of this process now. 




Figure 3: The cascade combines the coordinates of the four black 
circles with and 1 to form the five white circles. 



Cascades. We let i + 1 be the number of extreme dots and, 
for a reason that will become clear shortly, denote their coor- 
dinates with different indices as (uk, Wk+i), for < k < I, 
ordering them such that itrj < u\ < . . . < U(. By defini- 
tion of extreme, the corresponding second coordinates sat- 
isfy io i < u)2 < . . . < wt+i. Observe that uq is the globally 
minimum value and wi+\ is the globally maximum value. 
Let us now construct the reduced diagram, focusing on di- 
mension 0. We start the upward pass with the birth of the 
(— l)-dimensional homology class at 0, which dies at uq, 
so we have (0, ito) in Dgm_ 1 (g). The minima of all other 
components mark the births of 0-dimensional classes. They 
can be interpreted as gaps between components. We start 
the downward pass with the birth of a 0-dimensional class 
at 1, which dies at We+i, so we have (1, we+i) in Dgm (g). 
The maxima of all other components mark the deaths of 0- 
dimensional classes. If w is the second coordinate of a non- 
extreme component, then it marks the death of the gap that 
opened up when we passed the corresponding first coordi- 
nate, it, during the upward pass. Hence, (it, w) belongs to 
Dgm (g), and we note that this dot also belongs to the stan- 
dard diagram. However, if w = Wk is the second coordi- 
nate of a extreme dot, then it marks the death of the gap 
that has opened up when we passed the minimum of the next 
component with extreme dot during the upward pass. This 
minimum value is Uh, so (itfe, iffc) belongs to Dgm (g), for 
1 < k < I 
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In summary, we see that contiguous extreme dots ex- 
change one pair of coordinates, using on the left and 1 
on the right to complete the process. More precisely, the re- 
duced diagram is the same as the standard one, except for 
substituting the dots (JTvJ for the dots ([Tol l: 

(u ,wi), (i4i,to 2 ), • ■ • , (ue-i, wt), (ue,we+i), (16) 
(0, it ), (ui,tui), (u2,w 2 ), • ■ • , (u£, Wi), (1, w £+ i). (17) 

For later use, we call this substitution a cascade and write 
Dgm(g) = Dgm(g) c ; see Figure [3] Note that all dots be- 
long to the horizontal and vertical subdiagrams of dimension 
0, except for (l,we+i), which belongs to the relative sub- 
diagram of dimension 0, and (0, uq), which belongs to the 
ordinary subdiagram of dimension —1. We will often have 
uq = and W£+i = 1, in which case the first and the last 
dots in (fTTI i lie on the diagonal and can be ignored. 

Point calculus. As described in |[T] |2), the persistence di- 
agram can be harvested for a wealth of homological infor- 
mation. This includes the ranks of the homology groups of 
the sets X t and of the pairs (X, X*). To explain this, we write 




(M) (M) 

Figure 4: The persistence diagram of the function g — /|m, as de- 
fined in FigureQ] the circle dots also belong to the diagram of f\u, 
but the square dots do not. The rectangle on the left, Ct(g), repre- 
sents the homology of Mt, while the rectangle on the right, IZt (g), 
represents the homology of (M, M ); compare with the Betti num- 
bers computed for the example at the end of Section|2] 

£t(g) for the multiset of dots in the rectangle with lower cor- 
ner (t,t) on the base edge of Ord p (_g). Similarly, we write 
~R.f(g) for the multiset of dots if the lower corner, (t, t), lies 
on the base edge of Re\ p (g); see Figure|4] Often, we drop 
the dimension from the notation and write C, t {g) for the dis- 
joint union of the multisets £f (<?), over all p, and similar for 
1Zt(g). We can read the Betti numbers of X f by collecting 
the dots in C t (g), and we can read the relative Betti numbers 
of (X,X*) by collecting the dots in TZ t (g). Writing C p t (g) 
and IZt(g) for the corresponding multisets in the reduced 
diagrams, we note that the same relations hold between the 
reduced Betti numbers and the reduced persistence diagrams. 



Betti numbers and persistence diagrams. To motivate 
the first result of this paper, we rewrite ((T) to (O for reduced 
Betti numbers: 

&(V t ) = /yU,U 4 ), (18) 

for all < p < n and all regular values t. Observe also 
that V t = iff U - U* = § - §*, which implies that <0 
also holds for p = —1. In words, for every dimension p and 
every regular value t, the rectangles £P t (/|v) and 72f (/|u) 
have equally many dots. This 1-parameter family of rela- 
tions is satisfied if the reduced diagrams of /|v and /|u are 
reflections of each other. This is our first result, stated as the 
Land and Water Theorem in Section [4] We see an illustra- 
tion in Figure |2] which shows the reduced and non-reduced 
diagrams of the function / in Figure Q]restricted to the solid 
torus, U, and to the complementary solid torus, V. Remov- 
ing the black dots, we are left with the two reduced diagrams, 
which are indeed reflections of each other. 

To motivate our second result, we consider the relations 
( fT3l ) to dl~5b . They say that for every regular value t, the rect- 
angle £f (/|m) has the same number of dots as TZj (/|u) and 
^-t (/Id) together. Similarly, Vf t (/|m) has the same number 
of dots as L\ (/|u) an d T^-t(f\v) together. This 1-parameter 
family of relations is satisfied if the persistence diagram of 
/|m is the disjoint union of the diagram of /|u and of its re- 
flection. This is our second result, stated as the Euclidean 
Shore Theorem in Section [5] It is illustrated in Figure [4] 
which shows the persistence diagram of the function /|m in 
Figure Q] Comparing Figures [2] and [4] we see that the two 
circle dots in Figure [4] also belong to the diagram of /|u, 
whereas the two square dots are reflections of the circles. 

While ([]]) to © and ( fT3T > to <IT3T > motivate our two results, 
these relations are not sufficient to prove them. Indeed, ( fT3l ) 
to (fT~5b hold in general, but the Euclidean Shore Theorem re- 
quires that the minimum and maximum of the perfect Morse 
function belong to a common component of V. If this con- 
dition is violated, then we can have coordinate exchanges 
among the dots that contradict the Euclidean Shore Theorem 
without affecting the relations ( fT~3T > to ( fT~5T >. An example of 
this phenomenon can be seen in Figure |9j which will be dis- 
cussed in Section[5] 

4 Land and Water 

In this section, we present our first result, which extends 
Alexander duality from spaces to functions. Its proof needs a 
general algebraic result about contravariant filtrations, which 
we present first. 

Compatible pairings. A pairing between two finite-di- 
mensional vector spaces X and Y over the field F is a bi- 
linear map (, ) : X x Y — > F. The pairing is non-degenerate 
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if every ( £ X has at least one 77 G Y with (£, 77) = 1, and, 
symmetrically, every element 77 G Y has at least one ^ £ X 
with (£, 77) = 1. The term is justified by the fact that a non- 
degenerate pairing implies an isomorphism between the vec- 
tor spaces. To see this, we choose arbitrary bases of X and Y, 
writing n\ and ny for their sizes, and represent the pairing 
by its values between the basis elements. This gives an in- 
vertible 0-1 matrix with n\ rows and ny columns. Reducing 
the matrix to diagonal form gives the claimed isomorphism. 

Now consider two filtrations, each consisting of m + 1 
finite-dimensional vector spaces over F, and non-degenerate 
pairings connecting the filtrations contravariantly: 



X - 


■> X X -> . 


— > X m 


X 


X 


X 


■ m ^ 


- Y m _i . 


• ■ <- Y . 



It will be convenient to assume that the filtrations begin and 
end with the trivial vector space: Xo = X m = Yo = Y rn = 
0. Write g\ : X £ ->■ X., and hj : Yj ->■ Yj for the maps 
upstairs and downstairs. We call the pairings compatible 
with these maps if (£,77) = (£',77') f° r every £' G X^ and 
r? G Y m _j, where £ = gf(£') and 77' = h^.(r?). By what 
we said above, we have isomorphisms isOj : Xj — > Y m _j, 
for all < j < m, but we need more, namely isomor- 
phisms that satisfy a compatibility condition with the hori- 
zontal maps. Before defining what we mean, we prove a few 
properties of the compatible pairings. To begin, we split the 
vector spaces connected by gj : 

Xi = pimg; kerg|, (19) 
Xj = img; cokg;'. (20) 

and similarly for Y m _j and Y m -i', see Figure [5] Here, we 
have of course pim = X^/kerg^ and cokg^ = X^/im g 3 -. 
The restriction of the pairing to the kernel of g\ times the 




Figure 5: Splitting the spaces in the two anti-parallel filtrations. 



image of h™Z*- is well-defined because elements of the ker- 
nel and the image are also elements of the vector spaces. In 
contrast, an element of a quotient is an entire coset, whose 
members may have inconsistent values under the pairing. To 
define the restriction, we assume a basis that splits along the 
same lines, so that the coset in pim g\ has a unique represen- 
tation of the form £' + kerg^, in which the projection of £' 
to the kernel is zero. Finally, we define the restricted pair- 
ing using £' as a proxy for the coset it defines. Similarly, we 
define the restriction of the pairing to the cokernel of g\ . We 
have the following properties: 

(A) The pairing restricted to im g\ x ker h™I* is trivial, that 
is: all pairs map to 0. 

(B) We can find a basis of the cokernel such that the pairing 
restricted to cokg^ x pim k™Z) is trivial. 

To prove (A), let £ G im g\ and 770 G ker h™Z} both be non- 
zero. The image of 770 is 0, and a preimage, £', of £ exists. 
This implies (£, 770) = (£', 0) = 0, as claimed. 

To prove (B), let £0 G Xj and 77 £ Y m -j both be non- 
zero, with zero projection of £0 to im g\ and of 77 to ker h™I* ; , 
Recall that £0 and V act as proxies for the cosets £0 + im g^ 
in cokg^ and 77 + kerh™!} in pim h™Z*-. If (&,»?) = 0, 
then there is nothing to do, else we take steps to set the value 
of the pairing for the cosets to zero. Assume without loss 
of generality that £0 is an element of the chosen basis of Xj . 
By the non-degeneracy of the pairing X; x Y m _i, there exists 
£' G Xj with (£',77') = 1, where 7/ = h^I}(r?). By (A), 
we can assume that the projection of £' to the kerg^ is zero. 
By compatibility, we have (£,77) = 1, where £ = gf(£')- 
Now we have (£,o,v) = (ZiV) — 1> and £0 7^ £ because 
one belongs to the cokernel and the other to the image of 
g\. Replacing £0 by £1 = £0 — £ in the basis of Xj, we 
get (£1,77) = by bilinearity. We note that we still have 
a basis that splits as required. The operation achieves its 
goal by changing the restricted pairing without affecting the 
pairing. Repeating the operation for other basis elements of 
the cokernel, we eventually get a basis such that the pairing 
restricted to cokg^ x pim h™I* is trivial, as claimed. 

Compatible isomorphisms. We are now ready to specify 
the compatibility condition for the vertical isomorphisms. 
First, we require that they split along the same lines as the 
spaces. Specifically, we consider indices 0<i<j<k< 
m and the maps gj and gj, and we define 

X* = img; kerg;. (21) 

X 1 ;' = img; X k . (22) 

Xf = korg ! : Xf, (23) 

X= p = X,/(X^©X;©Xf); (24) 
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see Figure[6] For the case i + 1 = j = k — 1, these spaces 



because X^ k 



mg j andY^.©Y^ = kerhp:}. 
The lower-left four blocks are zero because X^ k © X^ p = 

cok gj and Y c £ l _ J © Y'm-j = P im n m-}- Repeating the argu- 
ment with reversed roles (swapping Xj with Y n -j), we con- 
clude that the remaining four non-diagonal blocks are also 
zero. Since the entire matrix is invertible, this implies that 
all four diagonal blocks are invertible. We finally get the 
isomorphisms by turning the diagonal blocks into unit ma- 
trices. All we need here are standard row and column oper- 
ations, which preserve the pairing while changing the basis 
for its description. The compatibility condition < f29b is there- 
fore but a reformulation of the compatibility condition for 
the pairings. This completes the argument for the existence 
of compatible isomorphisms. The reason for our interest is 
the following straightforward result. 

CONTRAVARIANT PE THEOREM. Let X to X m and Y 
to Y TO be two filtrations contravariantly connected by com- 
patible isomorphisms. Then there is a bijection between the 
classes upstairs and downstairs such that a class is born at 
Xi and dies entering Xj iff the corresponding class is born at 
Ym-j+i and dies entering Y m - i+ x. 

Proof. Recall that isoj : Xj — > Y m _j splits into four iso- 
morphisms, for each choice of triple indices i < j < k. The 
case i + l= j = k— 1 implies that the births and deaths 
upstairs correspond to the deaths and births downstairs. Con- 
dition (f29t implies that the births and deaths are matched up- 
stairs the same way as downstairs, only backward. ED 




Figure 6: Splitting the vector spaces into four. 



represent the classes that end at Xj, that extend before and 
after Xj, that exist only at Xj, and that begin at Xj, in this se- 
quence. Similarly, we consider the anti-parallel maps, h™~^. 
\Zj, and we define the subspaces Y^_j, Yj^ •, YJ^_j, 
downstairs. We call the vertical isomorphisms com- 



and ^ 

m-j 

patible with the horizontal maps if isOj 
into isomorphisms 



Xj ->■ V m -j splits 



ISO 



ISO 



ISO 



x;- 

XJ- 

xf 

Xf 



m-ji 
Y'P 

m-j i 
T m-ji 

. Y ik _- 



(25) 
(26) 
(27) 
(28) 



and the maps between the subspaces compose to the identity, 
for every < j < to; see Figure [5] More formally, this 
means that the maps satisfy 



<■ = ^(isor>™4(is Oj (0))), 



(29) 



for all £ G im g\ and < i < j < to. 

We now prove that the compatible pairings imply the ex- 
istence of compatible isomorphisms. The proof is construc- 
tive. To start, we choose bases of Xj and Y m -j that split four 
ways, as do the spaces. Representing the pairing between Xj 
and Y,„_j by the thus defined 0-1 matrix, we get a decom- 
position into 16 blocks; see Figure[7] We use (A) and (B) to 





Tit —j 





















xj 













xf 













xf 














Figure 7: The pairing between Xj and Y m _j splits into four non- 
degenerate restrictions. 

show that the bases can be chosen so that all non-diagonal 
blocks are zero. Indeed, the upper-right four blocks are zero 



We see that the compatibility condition d29l ) is for con- 
travariantly connected filtrations what commutativity is for 
covariantly connected filtrations. While the (covariant) 
Persistence Equivalence Theorem f3] page 159] implies 
equal persistence diagrams, the Contravariant P(ersistence) 
E(quivalence) Theorem implies reflected persistence dia- 
grams. 

Alexander duality for functions. After getting the ab- 
stract prerequisites settled, we are ready to state and prove 
our first result. Let / : S — > [0, 1] be a perfect Morse func- 
tion on the (n + l)-dimensional sphere, and § = UUYa 
decomposition into complementary subsets whose intersec- 
tion, M = U n V, is an n-manifold. Assuming /|m is tame, 
we have only finitely many homological critical values, in- 
cluding and 1. Letting these critical values be sq to s m , 
we interleave them with homological regular values to get 
So < to < Sx < . . . < t m _x < s m . We are interested in 
the filtrations defined by the restrictions of / to U and to V. 
To describe them, we write Of for the p-th reduced homol- 
ogy group of V ti , if < i < m, of U, if i = m, and of 
(U, U* 2m_i ), if m < i < 2m. In the same way, we define 
for < % < 2m with respect to V. Using Alexander duality 
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to form isomorphisms, we get two contravariantly connected 
filtrations: 



Uo 



V 



— > 



u 



v; 



V 



where p + q = n. We will prove shortly that the isomor- 
phisms implied by Alexander duality are compatible with the 
horizontal maps. Assuming this much, we get our first result, 
which we formulate using the superscript 'T' for the opera- 
tion that reflects a dot across the vertical axis of a persistence 
diagram and, at the same time, changes its dimension from p 
to q = n — p. 

Land and Water Theorem. Let U and V be comple- 
mentary subsets of § = S™+\ and let / : S [0, 1] be a 
perfect Morse function whose restriction to the n-manifold 



U n Vis tame. Then Dgm(/ 



Dgm(/| 



PROOF. By the Contravariant PE Theorem, we have a bisec- 
tion between the classes upstairs and downstairs that respects 
the matching of births and deaths. Specifically, if a class 
a is born at Of and dies entering 0^ then its correspond- 
ing class is born at V| TO _,- +1 and dies entering V| m _ i+1 . If 
1 < i < j < m, then the class a is represented by (si, Sj) 
in the p-th ordinary subdiagram of f\jj. The correspond- 
ing class is represented by (sj, Si) in the q-th relative sub- 
diagram of /|v- The reflection maps the first and second 
coordinate-axes of the ordinary subdiagram to the second 
and first coordinate-axes of the relative subdiagram. It fol- 
lows that it maps (s^ Sj) in the former to (sj 7 Si) in the latter 
subdiagram. Other cases are similar, and we conclude that 
Dgm p (/|u) and Dgm 9 (/|v) are reflections of each other. 
Writing this more succinctly gives the claimed relationship 
between the reduced persistence diagrams of / restricted to 
the two complementary subsets of the sphere. ED 



Alexander pairing. We fill the gap in the proof of the Land 
and Water Theorem by establishing compatible pairings be- 
tween the groups upstairs and downstairs. For the sake of 
simplicity, we restrict ourselves to the case F = Z2. While 
the argument for general fields is similar, it requires oriented 
simplices, which is a technical formalism we prefer to avoid. 
We begin with the pairing implicit in Lefschetz duality for 
manifolds with possibly non-empty boundary: 

(,) L : H p (U t ) x H, +1 (U tj 0U t ) ^Z 2 , 

which is defined by mapping two classes to the parity of the 
number of intersections between representing cycles. We get 
such a pairing for every regular value t, and these pairings are 
compatible with the horizontal maps induced by inclusion of 



the sublevel sets; see 0. Next, we reduce the dimension of 
the second factor using a mapping, 92, which we compose 
from four simpler mappings: 

H 9+1 (u t ,au t ) n H g+1 (S,§* UVt) 

<P 2 



^ H g (S*UV t ) 



H g (S'uV t ,S*) 



H,(V,V). 

The first mapping, ip\, is an isomorphism defined by exci- 
sion. Next, cp2, is the connecting homomorphism of the ex- 
act sequence of the pair (S, § f U Vt). Since the q-th reduced 
homology of § is trivial, for all q 7^ n + 1, 1^2 is an isomor- 
phism for all < q < n. It is a surjection for q = n. The 
third mapping, 933, is induced by inclusion. It occurs in the 
exact sequence of the pair (§ f U Vt, §*). For < t < 1, the 
reduced homology groups of S* are all trivial, implying that 
933 is an isomorphism. For t = 0, 993 is the trivial isomor- 
phism. Finally, tp 4 is again an isomorphism defined by exci- 
sion. In summary, we get ip : H g +i(Ut,5Uf) — > H 9 (V, V), 
which is an isomorphism for < q < n and a surjection 
for q = n. More specifically, each class (3 £ H ra (V, V*) has 
two preimages. Indeed, an n-cycle representing the corre- 
sponding class in H„(S* U Vt) partitions the components of 
U< into two subsets, and each subset generates an (n + 1)- 
dimensional relative homology class that maps to f3. 

Note that the roles of U and V can be interchanged to get 
a mapping from H q+1 (W t ,d¥ t ) to H,(U, U*). With this, we 
are ready to construct the Alexander pairing: 



(,) A :H p (U t )xH g (V,V*) 



So 



defined by (a,f3) A = (a,(3') L , where /3' is a preimage of 
(3 under (p. This fixes the pairing of Of and ^m-i f° r — 
i < m. Similarly, we define (,) A for H p (U,U*) x H 9 (V t ) 
which fixes the pairing for in < i < 2m. The pairing is 
clearly well-defined for < p < n, where cp is an isomor- 
phism. The remaining two cases are symmetric, and we con- 
sider p — and (3 G H„(V, V). As noted before, /3 has 
two preimages, f3' and f3", which are generated by comple- 
mentary subsets of Vt- Every a G Ho(Ut) is represented by 
an even number of points, which the complementary subsets 
partition into two even or two odd subsets. In either case, we 
have (a,f3') L = (a,/3") L , so the pairing is well-defined in 
all cases. Moreover, the compatibility of the Alexander pair- 
ings within the left and right halves follows from the com- 
patibility of the Lefschetz pairings. 

It remains to define the paring in the middle of the se- 
quence, for 0^ = H P (U) and V?, = H 9 (V). To this end, let 
a be a reduced p-cycle in U, f3 a reduced q-cycle in V, and 
/3' a (q + l)-chain whose boundary is /3. We define (a, (3) A 
by counting the intersections between a and /?'. Note that f3' 
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is the preimage of j3 under the previously defined mapping, 
if : H 9+ i(U t ,0U t ) -> Hq(V,V*), for any regular value t. 
Therefore, the pairing in the middle is compatible with the 
pairings in left half, again exploiting the compatibility of the 
Lefschetz pairing. Alternatively, we can define the pairing 
by taking a (p + l)-chain a' with boundary a and count- 
ing its intersections with j3. In this case, a' is the preimage 
of a under the map tp : H p+ i(V t , 9V t ) -> Hp(U,U*), and 
compatibility with the pairings in the right half follows. In- 
deed, both definitions are equivalent, as already observed by 
Lefschetz |8 j. We give a simple proof for the case F = Z2: 

Bridge Lemma. Let a and (3 be non-intersecting re- 
duced cycles on the (n + 1) -sphere whose dimensions add 
up to n. Let a 1 and j3' be chains whose boundaries are a and 
j3, respectively. Then (a', (i) L = (a, f3') L . 

Proof. We can assume that a and f3 do not intersect. The 
intersection of a' and f3' is a 1 -chain. Its boundary consists 
of an even number of points, and is the disjoint union of in- 
tersections of a' with j3 and of a with j3'. This implies that 
both types of intersections occur with the same parity. ED 



5 Shore 

This section presents our second result, which extends the 
combination of Alexander duality and Mayer- Vietoris se- 
quences from spaces to functions. 

Mayer- Vietoris sequence of nitrations. Assuming that 

/|m is tame, we write s < t < si < . . . < i m _i < s m 
for the interleaved sequence of homological critical and reg- 
ular values, as before. The main tool in this section is the 
diagram obtained by connecting the filtrations of / and its 
restrictions with Mayer- Vietoris sequences. We describe this 
using shorthand notation for the groups. Consistent with ear- 
lier notation, we write Sf for the p-th homology group of S ti , 
for < i < m, of S, for i = m, and of (S, §* 2m_i ), for 
m < i < 2m. Similarly, we write M P for the p-th homology 
group of M ti , of M, and of (M, M* 2m - ; ). Finally, we write 
D P for the direct sum of the p-th homology groups of and 
V tl ,ofUandV, and of (U, U' 2 —*) and (V,V* 2 "-*). Draw- 
ing the filtrations from left to right and the Mayer- Vietoris 
sequences from top to bottom, we get 



I 


I 


1 


s p+1 -> 


v CP+l V 

T J m T 


, CP+1 

... — f o 2m 


1 


I 


1 


M p -> ... 


Mf„ -> 




1 


1 




D p -> ... 


-> Df n -> 




1 


1 


1 


S p ... 


->■ S p -> 


• ■ • ^ °2m 


1 


1 


1 



All squares commute, which is particularly easy to see for the 
squares that connect groups of the same dimension, whose 
maps are all induced by inclusion. For 1 < p < n, all groups 
S P are trivial. By exactness of the Mayer- Vietoris sequences, 
this implies that the maps m p : M P — > D P are isomorphisms, 
for 1 < p < n — 1. The persistence diagram of /|m is 
therefore the disjoint union of the persistence diagrams of 
/|u and /|v- More generally, we claim: 

General Shore Theorem. Let n be a positive inte- 
ger, let U and V be complementary subsets of § = S n+1 , 
and let / : S — >• [0, 1] be a perfect Morse function whose 
restriction to the ?i-manifold M = U n V is tame. Then 



Dgm (/| M ) = [DgmoC/MuDgmot/hOf, 
Dgm p (/| M ) = Dgm p (/|u) U Dgm p (/| v ), 
Dgm„(/| M ) 



(30) 
(31) 



[Dgm (/|u) U Dgm (/| v )] Gi , (32) 

for 1 < p < n — 1, where C stands for applying the cascade 
and T stands for reflecting the diagram. 

We note that the assumption of n being positive is necessary 
since the formulas do not hold for n = 0. Concerning the 
proof of the theorem, note that ( Bit is clear, and that d32b 
follows from (f30b using the duality of persistence diagrams 
from [3 |. We will need to study the impact of the non-trivial 
groups to prove ( f30l >. 

Latitudinal manifolds. Call the minimum of / the south- 
pole and the maximum the north-pole of the sphere. Let 
M' be a component of M, and note that it is an n-manifold 
that decomposes § into two complementary subsets. If it 
separates the two poles, we refer to M' as a latitudinal n- 
manifold. Assuming neither pole lies on M, we order the lati- 
tudinal ?i-manifolds from south to north as Mi , M2, . . . , M.g. 
Letting Uk and Wk be the minimum and maximum values of 
/ restricted to Mfc, we get u\ < 112 < ■ ■ ■ < ui as well as 

Wi < U>2 < . ■ ■ < Wg. 

For each component M' of M, there are neighboring com- 
ponents U' of U and V of V defined such that D'(lV' = M'. 
A single component of U or V can be neighbor to an ar- 
bitrary number of n-manifolds, but not to more than two 
latitudinal n-manifolds. Specifically, there are components 
§0, §1, ...,§*: of U and V such that M fc = § fe _i n § fe for 
each 1 < k < £. We refer to the Sfe as latitudinal compo- 
nents. For example, in Figure [9] we have I = 2 latitudinal 
1-manifolds and £+1 = 3 latitudinal components. Setting 
Mo = and Wg+i = 1, the minimum and maximum values of 
/ restricted to are and Wk+i, for < k < I. Assum- 
ing Sfe belongs to U, it gives rise to zero or more dots in the 
ordinary subdiagram of Dgm (/|u) and to exactly one dot, 
(v,k, Wk+i), in the horizontal subdiagram. We say the dot in 
the horizontal subdiagram represents the 0-dimensional ho- 
mology class defined by E>k- Note that its coordinates are 
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indexed consistently with the notation used in the introduc- 
tion of the cascade. There is indeed a connection, namely the 
dots in the 0-th horizontal subdiagram representing latitudi- 
nal components are the extreme ones in the multiset: 

Extrema Lemma. The latitudinal components of U and 
V correspond bijectively to the extreme dots in the disjoint 
union of Dgm (/|u) and Dgm (/| v ). 

We omit the proof, which is not difficult. The statement in- 
cludes the case in which there is no latitudinal ?i-manifold 
so that §o is the only latitudinal component. It contains both 
poles and is therefore represented by (0, 1), which is the only 
extreme dot in the disjoint union of 0-th diagrams. 

Proof of d30l l. We compare Dgm (/|M) with the dis- 
joint union of Dgm (/|u) and Dgm (/|y), noting that their 
dots all belong to the ordinary and horizontal subdiagrams. 
Consider first a dot (a, b) in the ordinary subdiagram of 
Dgm (/|jyii)- It represents a component in the sublevel set 
that is born at a minimum x G M, with f(x) = a, and that 
dies at a saddle point y G M, with f(y) = b. Let M'CMbe 
the connected n-manifold that contains x and y, and recall 
that M'=lI'nV. Assume the neighboring component U' 
of U lies above x. In the sequence of sublevel sets of /|u, 
we see the birth of a component at f(x) = a and its death 
at f(y) = b. It follows that (a, b) is also a dot in the ordi- 
nary subdiagram of Dgm (/|u) and therefore of the disjoint 
union of the diagrams of f\\j and /|y. The argument can be 
reversed, which implies that the 0-th ordinary subdiagrams 
are the same. 

Consider second a dot (a, b) in the horizontal subdiagram 
of Dgm (/|ia). It represents a connected n-manifold M' C 
M, which splits § into two subsets. If M' is non-latitudinal, 
then one subset contains both poles and the other contains 
neither. The latter subset contains a (non-latitudinal) neigh- 
boring component, which is represented by (a, b) in the hor- 
izontal subdiagram of Dgm (/|o) or of Dgm (/|v)- Again, 
the argument can be reversed. If on the other hand, M' is lat- 
itudinal, then (a, b) = (ui,Wi), for some i, where we write 
iti, U2, ■ ■ ■ ,ut for the minimum values and w\, u>2, ■ ■ ■ , wi 
for the maximum values of the latitudinal ?i-manifolds, as 
before. More generally, we get the dots (itj, «;,•) in the hor- 
izontal subdiagram of Dgm (/|M), for 1 < i < i. Adding 
uq = and W( + i = 1, we get the dots (itj, lOj+i) in 
the horizontal subdiagram of Dgm (/|u) U Dgm (/|v), for 
1 < i < £. After adding (0, 0) and (1,1), which are both 
irrelevant, we get precisely the dots specified in ( fTTb and in 
dT6b . This implies that the two diagrams are related to each 
other by a cascade, which completes the proof of d30l >. 

As an alternative to the above geometric proof of ( f30l >. 
we could give an algebraic argument that considers the fil- 
trations of Meyer- Vietoris sequences; see Figure [S] which 



S 1 : : : : : 

M° : 

M- 1 i i i i i : 

OabcdldcbaO 

Figure 8: White areas represent trivial groups. The pattern of hav- 
ing every third row white is broken at the top and at the bottom. 
The black life-time intervals show homology classes that appear in 
the analysis of the 0-th and n-th persistence diagrams. 

draws the nitrations as rows in a matrix, shading the area of 
non-trivial groups. The obstacle to M" and being isomor- 
phic, and thus the reason for the cascade, are the non-trivial 
groups S? in the left half of the filtration. Because each non- 
trivial S° has exactly one generator, must have one more 
generator than M". Indeed, by comparing the groups, before 
and after the cascade, we see that in the left half, each has 
one more generator than M°. 




Figure 9: Left: the height function on the 2-sphere decomposed 
into an annulus and a pair of disks. Right: the reduced persistence 
diagram of / restricted to the annulus (shaded circular dots), to the 
pair of disks (white circular dots), and to the two circles (squares). 

Euclidean case. Considering that (|T3T > to ( fT3T > give elegant 
relations between the Betti numbers of the shore and the 
land, it is perhaps surprising that we need cascades to for- 
mulate a similar result for persistence diagrams. Indeed, ( TO! ) 
to ( fT~5T > suggest that the persistence diagrams of /|m be the 
disjoint union of the persistence diagram of /|u and of its 
reflection. The example in Figure [9] shows that this simple 
relation does not hold in general. Indeed, we have the dots 
(a, c), (b, d), (c, a), (d,b) in the standard diagram of /|m, 
while the standard diagram of /|u contains (a, d) and (c, b). 
This clearly violates the suggested relation. However, there 
is a natural setting in which the relation does hold, which we 
now describe. 

Euclidean Shore Theorem. Let n be a positive inte- 
ger, let A be a compact set whose boundary is an n-manifold 
in R™ +1 , and suppose that e : W l+1 — > K has no homolog- 
ical critical values and its restriction to dA is tame. Then 
Dgm(e| 9A ) = Dgm(e| A ) U Dgm(e|A) T . 
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Proof. We can extend e to a perfect Morse function / : 
§ n+i _^ K an d A to a subset U of § = § n+1 such that the 
persistence diagrams of e restricted to A and to <9A are the 
same as those of / restricted to U and M = <9U. It thus 
suffices to show that the persistence diagram of /|m is the 
disjoint union of the diagram of f\\j and of its reflection. For 
< p < n, this follows from (13 It . from 

Dgm p (/| v ) = Dgm„_ p (/|u) T , 

as stated in the Land and Water Theorem, and the fact that the 
reduced diagrams are equal to the non-reduced ones. For p = 
0, we start with ( l30t and note that there is only one extreme 
component in Dgm (/|u) U Dgm (/|v), namely the one of 
V that contains both the minimum and the maximum of /. It 
follows that the cascade leaves Dgm (/|u) unchanged while 
it turns Dgm (/|y) into Dgm (/|y). Using the Land and 
Water Theorem again, we obtain: 

Dgm (/| M ) = [Dgm (/|u) U Dgm (/| v )] c 

= Dgm (/|u) U Dgm (/| v ) 

= Dgm (/|u) U Dgm„(/|u) T 

= Dgm (/| u )UDgm„(/| u ) T . 

Finally, for p = n, we exploit Poincare duality for manifolds, 
which implies Dgm„(/ |m) = Dgm (/| M ) T ; see 0. □ 

The 4-dimensional version of the theorem, for n = 3, 
brings us back full circle to the motivation for this work, 
namely the computation of the persistence diagram of the 
space-time shape formed by a moving collection of biologi- 
cal cells J6). Modeling space-time as M 4 , the shape is com- 
pact, and we consider the time function restricted to that 
shape. This data satisfies the assumptions of the Euclidean 
Shore Theorem, so we can infer the persistence diagram of 
the function on the boundary from the diagram of the func- 
tion on the solid 4-dimensional shape. 

6 Discussion 

The main contributions of this paper are two extensions of 
Alexander duality from spaces to functions. The first exten- 
sion is direct and relates the persistence diagrams of a perfect 
Morse function restricted to two complementary subsets of 
the {n + 1) -dimensional sphere with each other. The second 
extension relates the persistence diagram of the function re- 
stricted to the intersection of the complementary subsets with 
the diagram of the function restricted to one subset. A key 
tool in its proof is the filtration of Mayer- Vietoris sequences 
(or the Mayer- Vietoris sequence of nitrations). This suggests 
the study of more general filtrations of exact sequences. Be- 
sides the hope to develop a general purpose device that per- 
vades persistent homology in the same way exact sequences 



pervade homology, we motivate the study with a few con- 
crete questions. 

• Is it possible to relate the extent to which a function / on 
§™ +1 is not perfect Morse with the severity with which 
this function violates our two extensions of Alexander 
duality? 

• Does every Morse function / on § = S n+1 have a de- 
composition into complementary subsets, § = U U V, 
such that up to small modifications, the disjoint union 
of the diagrams of / and /|m is equal to the disjoint 
union of the diagrams of f\u and /|y? 

• Can these or similar relations be developed into a 
divide-and-conquer algorithm for computing the persis- 
tence diagram of / on a sphere? 

It would furthermore be interesting to generalize the results 
of this paper, as well as the above questions, to functions on 
(n + l)-manifolds other than the sphere. 
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